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ABSTRACT 

We propose a general method for mapping the extinction in dense molec- 
ular clouds using 2MASS near-infrared data. The technique is based on the 
simultaneous utilization of star counts and colors. These two techniques provide 
independent estimations of the extinction and each method reacts differently to 
foreground star contamination and to star clustering. We take advantage of both 
methods to build a large scale extinction map (2.5° x 2.5°) of the North America- 
Pelican nebulae complex. With star counts and H — Kg color analysis the 
visual extinction is mapped up to 35 mag. Regions with visual extinction greater 
than 20 mag account for less than 3% of the total mass of the cloud. Color is gen- 
erally a better estimator for the extinction than star counts. Nine star clusters 
are identified in the area, seven of which were previously unknown. 



Subject headings: ISM: clouds — dust, extinction — ISM: individual (North 
America Nebula, Pelican Nebula) — infrared: ISM — stars: formation 
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1. Introduction 

The knowledge of the extinction in molecular clouds is essential for a wide range of 
purposes. The extinction traces the dust distribution and thus provides an estimate of 
the column density of hydrogen that can be compared to the far-infrared or submillimeter 
emission from dust and to molecular emission lines. Fundamental properties are derived 
from these comparisons. Stepnik et al. (2001) have found evidence of dust evolution in a 
translucent cloud in Taurus using extinction from J-2MASS star counts and submillimeter 
data from the balloon-borne experiment PRONAOS. Kramer et al. (1999) have compared 
the C^''0 line with extinction derived from H — K color and have found a depletion of CO 
for Ay > 10 which indicates the interaction between dust an gas in high density regions. 
The dust distribution also gives information about the fragmentation processes which may be 
responsible for the shape of the initial mass function (Padoan et al. 1997b) through turbulent 
motions in the interstellar medium (Padoan et al. 1997a). All these results were obtained 
with near-infrared data which are better adapted than optical wavelengths for extinction 
studies. 2MASS provides accurate (~ 3%) J, H and Kg photometry for the whole sky. The 
Kg band is about 10 times less absorbed than the V band making possible to analyze the 
dust distribution for visual extinctions as large as 40 mag. Since the dust both absorbs 
and reddens the starlight, the analysis of the star density and star color variations give 
two independent estimates of the extinction. In this paper, we take advantage of these two 
techniques to measure the extinction and understand better the limitations of each method. 

The first section of this paper describes the methods for mapping the dust distribution 
based on Kg star counts and H — Kg reddening. The comparison of the two maps leads to 
the identification of foreground stars and stellar clusters. In a second section, we apply our 
technique to the North America-Pelican complex which is located in the galactic plane at an 
intermediate distance (~ 580 pc). By analyzing the statistical and systematic uncertainties 
in the star count and in the color extinction maps we merge these two maps into a more 
rehable representation of the cloud. We also propose the identification of new stellar clusters 
in this complex field. In the appendix, we discuss the comparison between 2MASS data and 
a synthetic model for stellar populations. 

2. Extinction mapping 
2.1. Star counts 

Variation of the star density across the sky was the first observable indication of ex- 
tinction, i.e. of presence of dust in the Galaxy. Star counts are therefore the oldest method 
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of mapping the interstellar cloud extinction (Wolf 1923). The method consists of counting 
stars in magnitude intervals in each cell of a regular grid in an obscured field and comparing 
this number with the counts in a reference field, supposedly free of absorption. To improve 
the spatial resolution, Bok (1956) made counts up to the limiting magnitude rather than 
in magnitude intervals. Manual techniques were used on Schmidt plates for many decades 
(Dickman 1978; Cernicharo & Bachiller 1984; Mattila 1986; Gregorio Hetem et al. 1988; 
Andreazza & Vilas-Boas 1996). More recently Cambresy et al. (1997) proposed a variant of 
the classical method by replacing the regular grid by an adaptive grid. They used infrared 
data from the DENIS survey (Epchtein 1997) and an automatic counting technique that is 
possible with digital data. The adaptive grid optimizes the angular resolution with respect 
to the local star density. This is useful when there are significant fluctuations of density due 
to large extinction gradients or the presence of star clusters. 

In its simplest form, star count method assumes that all stars are behind the cloud and 
that the stellar population is homogeneous for the whole field. Moreover, the star density is 
supposed to be uniform within a cell. Cloud structures on angular scales smaller than the 
resolution of unit cell can lead to underestimates of the extinction (Rossano 1980; Thoraval 
et al. 1997). From star counts in the Ks band, the visual {V) extinction is derived from star 
density (D) as follows: 

Ak. = -log%^ = Z + ilogl (1) 
a D a D 

Ay = ^x^K. (2) 

where a is the slope of the Ks luminosity function and Dj-cf the density in a reference field. 
Since the reference density depends on the position (mainly on the galactic latitude), we use 
a zero point function, Z, which can be approximated by a constant value only for small areas 
(typically less than 1 deg^). For larger areas, Z is estimated using a fit of the uncalibrated 
Av versus galactic position for unobscured regions close to the cloud. 



2.2. Reddening 

Dust that obscures also reddens starlight. Measuring the reddening of stars behind 
a cloud gives the extinction along the line of sight. The spatial resolution of the derived 
extinction is very high and corresponds to the apparent size of the stellar disk, but the mea- 
surements are strongly undersampled. The reddening is the difference between the observed 
star color and the intrinsic color corresponding to the star spectral type. Frerking et al. 
(1982) used J — K color to map the extinction toward fields in p Ophiuchus and Taurus. 
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Jones et al. (1984) used the same color to map the extinction in a Bok globule toward the 
Southern Coalsack. Until recently, the necessity of knowing the spectral type precluded a 
wide use of this method. However, with the advent of large digital surveys like 2MASS, 
statistical studies become possible where the knowledge of the spectral type is replaced by 
the measure of the mean color in an unobscured field. Using 2MASS prototype data, Beich- 
man & Jarrett (1994) derived the extinction toward the Taurus cloud by comparing J — H 
and H — Kg colors to a galactic model of stellar populations. Lada et al. (1994) proposed 
a method based on the H — K averaged color in a regular grid to map highly extinguished 
clouds with visual extinction up to 30-40 mag (Alves et al. 1998; Lada et al. 1999). Using 
DENIS data Schultheis et al. (1999) mapped the extinction toward the galactic center with 
a sophisticated method based on the isochrone position in a J — Kg/Kg color-magnitude 
diagram. 

As for star counts, a homogeneous population is assumed for the reddening analysis. 
Stars must lie behind the cloud and the reddening is supposed to be uniform within a cell. 
For example, the study of Schultheis et al. (1999) in the direction of the galactic center was 
restricted to RGB/AGB stars; they selected a specific stellar population at a large distance 
(behind the galactic center). 

We have chosen to use the H — Kg color because it presents a small intrinsic dispersion 
over typical stellar spectral types. Typically, we have < H — Kg < 0.37 (BcsscII & Brett 
1988) and a single peak in the color distribution. The J — Kg color has separate peaks for 
dwarf and giant stars which makes it difficult to convert a mean color to an extinction. 

We propose here an adaptive variant of the Lada's method. The size of our cells change 
with the local stellar density and is varied to contain always the same number of stars 
detected in both H and Kg 2MASS bands. The motivation for using adaptive cells is the 
same as that mentioned for star counts: to optimize the angular resolution for large variations 
of extinction and increase the resolution around stellar clusters. Visual extinction is derived 
as follows: 

Er-Ks = {H - -ft's)obs - {H - ii's)int (3) 

where {H — Kg)ohs is the observed median color in the cell and (H — Kg)[^it is the observed 
median color in the reference field and represents the intrinsic median color for unreddened 
stars. This differs from the Lada's method in which mean color is used rather than median 
color. Our choice is driven by the foreground star correction that is discussed in section 
2.3.2. Color is less sensitive to galactic latitude than star counts and a constant characteristic 
value (e.g. zero point of extinction) is correct for several square degrees. On larger scales. 
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the dependence with the galactic latitude (and longitude) must be taken into account with 
a zero point which is a function of direction. 

2.3. Simultaneous use of the two methods 
2.3.1. Calibration 

To use both star counts and reddening to study a cloud it is necessary to ensure that 
the calibrations are consistent. First, the spatial resolution must be identical, i.e. same cells 
must be used. Practically, the cell size is adapted to contain a number of stars, N, detected 
in both H and Ks- The Ks star density is measured in the cell (number of Kg stars divided 
by the cell surface area) and the H — Kg color of the N stars is represented by their median 
value. 

Star counts and colors are converted into extinction using equations (l)-(4). Since the 
zero point of the color-based extinction map is more reliable than that of the star counts- 
based map due to a weaker dependence on galactic latitude, the star counts map is adjusted 
to the zero point of color map. The larger the area free from absorption used, the better is 
the adjustment. 

For a small field (~ 1 deg^), the difference between the maps is simply an offset. For a large 
cloud (< 100 deg^), a linear fit of the difference between the two maps versus the galactic 
latitude is required. For still larger sizes, variation with the galactic longitude may also occur 
and the galactic latitude fit may no longer be linear. 

2.3.2. Expected differences 

Once the extinction maps are constructed with a consistent calibration, significant dif- 
ferences are still expected because each method has its own behavior when the a priori 
assumptions are not valid. 

Effect of foreground stars 

Star counts. Foreground stars are not affected by the cloud extinction and their number 
density, , does not vary; whereas the background star density, N^, is affected by the 
extinction as described in equation (1): — Nq x 10"""^^^', where a is the Kg luminosity 
function slope and the subscript refers to zero-extinction stellar density. The observed 
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extinction can be written: 

= ^^og^T^-^ (5) 



1 NU 



1 Mb I AT f 

= -log ^ (6) 

Let X be the fraction of background stars, X — Nq/{Nq + N^), we have: 

^^bs ^ -llog(^X10-'^^S:' + (l-X)) (7) 

Foreground stars reduce the density contrast between low and high extinction regions. For a 
given X, the measured extinction cannot exceed — 1/a x log(l —X), even for an infinite real 
extinction. For example, with X — 0.9 the observed visual extinction cannot exceed 26 mag 
from Ks star counts (10.5 mag from J star counts and 3 mag from V star counts). 



Reddening. Since the color in a cell is represented by the median color of the stars it 
contains, the contamination by foreground stars makes the derived extinction drop to zero 
when more than half of the stars are in front of the cloud (i.e. when the median color 

corresponds to the foreground star color). This happens at the wavelength A for Ax — 
— l/alog[(l — X)/X]. In the H band, which is the most affected by extinction, the turn over 
occurs at Ay ~ 16 mag for X = 0.9. 



Foreground star suppression. In Figure 1, the observed extinction is plotted versus the 
true extinction for a fraction of background stars X. The star count plots are obtained using 
equation (7). The color excess plot is obtained assuming a color distribution characterized by 
H — Kg = 0.17 i: 0.1 mag that corresponds to the vahic obtained with the stellar population 
model described in the appendix, for a direction close to the galactic plane (/ = 85°, b = —1°). 
As expected. Kg star counts give better results for high extinction than J star counts because 
of the relative sensitivity to the extinction {Aj/Ak^ ~ 2.5). The color excess technique is 
definitely more robust to foreground star contamination than star counts until the number 
of background stars drops to equal the number of foreground stars. 

Consequently, for high extinction. Ay ~ 20 mag and X — 0.9, Kg star counts give 
a much higher extinction (f=i 15 mag) than does H — Kg reddening (almost mag). The 
star count extinction map would show an extinction peak where the reddening map would 
show a hole in the cloud. Using the two maps allows one to show that the extinction is so 
high that there are more foreground stars than background stars in a particular direction. 
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The foreground star density can be estimated from these directions using the J — Kg color, 
which is the most sensitive to extinction. Once the foreground star density is known, we can 
remove these objects for each direction assuming that they are the bluer (i.e. lower J — X,) 
population. This hypothesis is verified only for high extinction regions, where the reddening 
is larger than the total uncertainties in the two maps. For low extinction, removing the 
bluer stars introduces a bias that will be discussed later (§ 3.3). When foreground stars are 
removed from the catalog, the star counts and the reddening analysis must be iterated to 
build reliable extinction maps for which the foreground star contamination is corrected. 

Young stellar clusters 

Molecular clouds are star forming regions and young embedded star clusters are to be ex- 
pected in such areas. A cluster produces an artifact in the star count extinction map: the 
excess density is interpreted as a hole in a cloud. On the other hand, such clusters generally 
contain stars with an intrinsic infrared excess due to an envelope or an accretion disk (es- 
pecially for low-mass stars) which leads to an extinction peak m a. H — Kg color extinction 
map. 

This is the opposite situation that was encountered for the contamination by foreground 
stars: the extinction appears larger in the color map than in the star count map. Directions 
for which the star density is significantly too large compared to the color suggest the presence 
of star clusters. A similar approach is developed in Carpenter (2000) who identifies young 
clusters by comparing 2MASS Kg star counts to stellar densities predicted by a model toward 
Perseus, Orion and Monoceros R2 molecular clouds and using CO line to correct for the 
extinction. 



3. Application to the North America-Pelican nebulae 
3.1. 2MASS data 

The 2MASS survey provides near-infrared photometry in the J, H and Kg filters. At 
the time of this paper data covering, 48% of the sky has been released through the 2"^^^ 
Incremental Data Release (Cutri ct al. 2000) and 100% of the sky has been observed. Since 
data analyzed here arc not in the 2"^^^ release, wc have extracted the catalog directly from the 
point source working database which contains sources extracted in each Tile (a 2MASS unit 
of 8.5' X 6° obtained by the co-addition of 274 images). As distinguished from the 2"^* release 
Point Source Catalogue, the working database contains multiple entries for the same star in 



-8- 



Tile overlap regions, low signal-to-noise sources, and artifacts flagged during the processing. 
We have extracted from the point source working database sources that satisfy the following 
criteria which are similar to the criteria used to generate the point source catalogue: 

1. photometric uncertainty a < 0.25 

2. signal-to-noise greater than 7 in at least one band 

3. no filter glint flag (ghost sources that occur at specific geometric locations relative to 
bright stars) 

4. no diffraction spike flag 

5. persistence probability due an echo of a previously observed bright star < 0.5 in the 3 
bands simultaneously 

In this flltering, we kept sources bright enough to be accepted as a real source but weakly 
contaminated by a diffraction spike, and sources with high PSF fltting (generally due to 
confusion between two stars). We give preference to completeness rather than reliability 
because our analysis is statistical and involves star counts. 

Redundant source entries in the Tile overlaps were removed by flltering out sources 
within 3" of bright stars (brighter than 8.5 mag in one band) and within 1.5" of fainter stars. 
The flnal catalog contains about 1.1 x 10^ objects. 

Based on an estimate of completeness from the magnitude distribution. Figure 2, we kept 
only stars brighter than Kg < 15.0 mag for the star count analysis and stars brighter than 
H < 15.5 mag and Ks < 15.0 mag for the color analysis. Analysis of repeated observations 
of calibration fields observed ~ 100 times indicates that the integral 2MASS completeness 
for these magnitudes is ~ 99% (Cutri et al. 2000). The density and color distribution of the 
selected stars does not show any residual persistence of the 2MASS observing Tiles. This 
validates the catalog eoctraction from the database and the adopted limiting magnitudes. 
An incorrect treatment of multiple entries in the working database would have led to the 
detection of the Tile overlaps in the density map while an overestimation of the limiting 
magnitudes would have introduced a star density level difference between Tiles, especially 
for those observed at different dates, under different airmass and/or atmospheric conditions. 
To derive extinction from Kg star counts using equation (1) the slope of the luminosity 
function, a, is needed. Its estimate from Figure 2 is a — 0.34. 
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3.2. Cloud properties 

The cloud studied in this paper hes between two HII regions, the North America nebula 
(NGC 7000) and the Pelican nebula (IC 5070). This complex is also named W80 (Westerhout 
1958) and the dark cloud itself is known as L935 (Lynds 1962). In the following, we use 
NAN to refer to the whole complex. Visual extinction has been deduced from Ha surface 
brightness by Goudis & White (1979) who found a maximum visual extinction of about 
8 mag. Bally & Scoville (1980) performed ^^CO observations of the cloud and proposed 
an evolutionary scenario for the molecular cloud. They concluded that the cloud is highly 
fragmented and contains an expanding molecular shell. Wendker et al. (1983) proposed a 
model to described the whole complex based on radio continuum observations. They gave 
a distance estimate of 500 pc. This distance has been confirmed by Straizys et al. (1993) 
who obtained d = 580 pc using photoelectric photometry of 564 stars in the direction of 
the complex. The total mass of the cloud has been estimated from molecular observations 
(i^CO, ^=^C0 and H2CO) to be 5 x 10^ (Feldt & Wendker 1993). 

Herbig (1958) discovered 68 stars having the Ha line toward the complex, suggesting 
the existence of a population of Be (intermediate mass) and T Tauri (low mass) stars. Delia 
Prugna et al. (1984) reobserved these sources to provide more accurate astrometry. Welin 
(1973) and Marcy (1980) discovered additional Ha emission stars in the NAN. The total 
number of previously known young stars in the square area defined in the Figure 3 is 119. 

To clearly delimit the dark cloud, we performed the star counts and the color analysis on 
a large area around the cloud: 5.6° diameter centered on aj2ooo — 10''56™00^, 5j2ooo — 44°43™43' 
However, the extinction study is restricted to a 2.5° x 2.5° area centered on aj2ooo — 20*^54^05^, 
(5j2ooo = 44°13™57" {I = 84.9, b = -0.41). The optical image and J - map of the large 
region are presented in Figure 3. J — Kg color is more sensitive to low extinction than 
H — Kg and shows the general shape of the dark cloud. The reference field to set the zero 
point of extinction has been chosen in the area of lower J — Kg value. It contains 20, 000 
stars detected in both H and K^. 



3.3. Results 

In the following, we assume the Rieke & Lebofsky (1985) extinction law rather than the 
most widely used Cardelli et al. (1989) law because the former is in better agreement with 
the reddening vector slope measured in the 2MASS color-color diagram {H — Kg/ J — H) 
obtained for the NAN (Figure 4). We have measured a slope of 1.61 ± 0.32 and Rieke & 
Lebofsky (1985) and Cardelli et al. (1989) obtained 1.70 and 1.21, respectively. We will 
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therefore use the following values: AjjAy = 0.282, AnjAy = 0.175 and AkJ Ay = 0.112. 

3.3.1. Extinction 
Reddening method 

The extinction map of the NAN is obtained using H — Kg color with the method described 
in § 2.2. The average color in the reference field provides the absolute calibration of the 
extinction (zero point). The average color in the reference region is {H — Ks)hit = 0.16, with 
a dispersion (1 a) of 0.12 mag which corresponds to Ay = 1.9 mag. The dispersion represents 
the sum of the photometric uncertainties, the intrinsic dispersion due to the average over 
all spectral types and a possible contamination by a residual absorption is the field. This 
dispersion is consistent with the one obtained using the galactic model presented in the 
appendix (0.1 mag), suggesting that there is no significant extinction contamination in the 
reference field. To estimate the systematic error resulting from a possible residual absorption 
one can use the hydrogen column density. Molecular hydrogen column density can be derived 
from the CO emission. For this field W^(CO) = 0.03 K km s~^ (Dame et al. 2001) which 
corresponds to A^(H2) = 5.4 x 10^^ cm~^. The atomic hydrogen column density is obtained 
from the 21 cm emission line and is iV(Hl) = 3.7 x 10^^ cm~^ (Leiden/Dwingeloo survey. 
Burton & Hartmann 1994). This value represents the total hydrogen column integrated along 
the line of sight. A velocity analysis of the neutral gas allows one to separate the different 
galactic components on the light of sight (Figure 5). Following the Fountain et al. (1983) 
results based on a Ha study of the NAN, the LSR velocity of the nebula is —7.1 ±5.5 km s~^ 
with a FWHM of 28.6 ± 0.6 km s~^. The fraction of Hi emission coming from this velocity 
range represents 40 ± 7% of the total. This value is an upper limit since not all the gas in 
this velocity range is associated with the complex. For [A^(Hl) + N{}l2)]/Ay = 1.87 x 10^^ 
(Savage & Mathis 1979) the upper limit for the visual extinction in the reference field is 
about 0.78 ± 0.14 mag (only 0.4% is associated with the molecular gas). 

Star count method 

The extinction is derived from Kg counts as described in § 2.1. Once the color map has been 
calibrated (previous paragraph), the star count extinction map is adjusted using a linear 
relation for the zero point versus the galactic latitude and longitude. For the NAN area, 
the stellar density is found to increase by ~ 25% per degree toward the galactic plane. The 
cell size is allowed to vary to achieve the desired number of stars, 10 in this work. Since the 
H — Kg intrinsic color dispersion is low, this small number gives good results for the color 
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map {a{Av) — 1.9/VTO — 0.6 mag). However, this small number introduces a significant 
Poisson noise in the star count map in which the standard deviation of the extinction is 
a = 3.5 mag. For Ay ^ 20 mag, the uncertainty is actually reduced for Kg star counts 
because the cell size is constrained by the H band which is more sensitive to extinction. 
Cells in high extinction regions contain 10 stars simultaneously detected in H and Kg and 
more than 10 additional stars detected only in Kg. The statistical uncertainty in the Kg star 
count map is reduced to cr ~ 2.5 mag when Ay ^ 20 mag. 

Foreground star density estimation 

Figure 6 shows the extinction from star counts and from reddening in the most obscured 
portion of the field. The effect of the foreground star contamination is obvious: the extinction 
derived from star counts reaches 30 mag where the extinction derived from reddening drops 
suddenly to mag. In areas where Ay(starcounts) — Ay{color) > 20 (i.e. 5.5o") we find 
57 stars. Among them, we identify 15 stars which do not show any J — Kg excess (the 
expected excess when Ay = 20 is Ej_k^ = 3.4). The normalized density for this foreground 
star population corresponds to 1900 ± 500 stars deg~^ and represents ~ 6% of the total star 
number for an unobscured field. 

Foreground star correction 

Unfortunately, setting a color threshold to separate the two star populations systematically 
yields a foreground star distribution that is strongly correlated with the dust distribution. 
The degeneracy between background blue stars and foreground red stars prevents the proper 
separation of the two populations. To correct for the foreground star contamination, we re- 
move the 30 bluer (i.e. lower J — Kg) stars per 7.5' x 7.5' box in the whole field, corresponding 
to the foreground star density measured toward the densest part of the cloud. In high ex- 
tinction regions this procedure removes foreground stars, but for visual extinction lower than 
~ 15 mag the degeneracy prevents a correct separation of the two populations and a bias is 
introduced in the extinction maps. The resulting foreground star distribution is presented in 
Figure 7. The iteration of the star counts and the reddening analysis with the new version of 
the catalog provides the foreground star corrected extinction maps (Figure 8). The standard 
deviation of the difference between these maps and maps obtained before the foreground star 
correction are 0.1 and 0.9 mag for the color and the star count extinction map, respectively. 
The main difference comes from the cloud core which corresponds to a high extinction for 
both maps: visual extinction is larger than 25 mag in the color map where it was mag 
before the correction. 
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To estimate the error in the foreground star correction we compare the distribution of 
these stars with the dust distribution. The difference between the mean extinction in the 
whole field and the mean value measured at the position of each supposedly foreground star 
should give zero for uncorrelated distributions. We find AAy = —0.26 mag. Simulations 
for an uniformly distributed star field yield a standard deviation of a = 0.03. The spa- 
tial distribution of the selected stars is still correlated with the extinction and is therefore 
not consistent with a population of only foreground stars. The fraction of foreground stars 
in a 7.5' x 7.5' box reaches 0.5 for the darkest region and is less than 0.1 for Ay < 15. 
Consequently, in regions for which the stellar color degeneracy exists, fewer than 10% of 
the stars are affected by the selection bias. For higher values, the foreground star iden- 
tification is straightforward (i.e. not statistical). We conclude that the contamination by 
residual background stars in the sample has negligible consequences on the extinction map- 
ping compared to the intrinsic dispersion of these maps: the standard deviation between the 
corrected and the not-corrected maps arc 0.1 and 0.9 mag whereas the intrinsic uncertainty 
for the extinction maps are a = 0.6 and a = 3.5 mag for the color and the star count map, 
respectively. 



Merging the two maps 

The uncertainty is lower from the color analysis than for star counts for Ay < 15. For 

higher values, it is not clear whether star counts or color give the best answer. The color 
method is more constraining because of the necessity to detect stars in both H and Ks 
bands whereas only detections in Kg are required to build the star count map. Since stars 
arc detected only in the longer wavelength in the most extinguished areas, more data are 
available for the star counts than for the color analysis. Statistical uncertainties are lower 
in the star count map when the number of Kg stars is greater than 300 in a cell (that by 
construction always contains 10 stars detected in both H and Kg). Using the extinction law, 
the mean color H — Kg — 0.16±0.01 and the limiting magnitude difference for each band, 
^iim _ j^hra _ g g mag, it would corrcspond to Ay ~ 70 mag. Therefore, the color analysis 
technique has a lower statistical uncertainty everywhere in our field. 

The resolution of the extinction maps is adapted to the local star density. It is 1.2' 
for Ay ^ Q mag, 1.9' for Ay = 10 mag, 3.7' for Ay = 20 mag and 7' for Ay > 30 mag. 

Unfortunately, extinction is probably no longer homogeneous for resolutions coarser than 
0.5 pc (3' for d = 580 pc) and both methods underestimate the extinction in a non-linear 
way. This effect can be simulated using an uniform spatial distribution for the stars and 
a the luminosity functions observed with 2MASS toward the NAN direction (see figure 2). 
One can then redden this simulated stellar field with any known extinction profile, assuming 
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all stars are background and applying the 2MASS completeness cuts. The simulation of a 

Gaussian peak of 30 mag of visual extinction with a 50" FWHM observed at 100" resolution 
gives Ay{H — Kg) = 4.6 mag and Ay{Ks counts) = 8.0 mag. A 90" Gaussian peak observed 
at a 100" resolution gives Ay{H — Kg) = 16.5 mag and Ay{Ks counts) = 18.8 mag. When 
the extinction is not homogeneous over a cell, star counts lead to a higher value than red- 
dening. But both methods underestimate the peak value. Thus, at the highest extinctions, 
using these techniques at the 2MASS sensitivity underestimates the extinction. Deeper ob- 
servations are needed such as those presented in Alves et al. (2001) with hmiting magnitudes 
of about 20 mag in both H and Kg. 

Our analysis concludes that (1) the color method is more robust than stars counts to 
the foreground star contamination. This contamination can be corrected by comparing color 
and count maps. (2) The star count extinction map has larger statistical uncertainty than 
the color extinction map. (3) The star count analysis is more robust than the color analysis 
to extinction inhomogeneities within a cell, i.e. the systematic uncertainty is smaller in the 
star count map for high extinction. Because of these different characteristics, we decided to 
combine the two extinction maps using reddening for Ay < 15, star counts for Ay > 25 and 
a linear combination of the two maps for the intermediate extinction range: 

Ay ^xx Ay(Ks counts) + (l-x) X Ay(H - Ks) (8) 

where: 

X — for Ay{KsCOunts) < 15 

X = [Av{K, counts) - 15] /lO for 15 < counts) < 25 

X — 1 for Ay (Kg counts) > 25 

The final combined extinction map of the NAN is presented in Figure 9. The maximum 
visual extinction is Ay — 35 mag whereas previous works at optical wavelengths estimate 
the extinction to less than 10 mag (Goudis & White 1979; Straizys et al. 1993). The 2MASS 
near-infrared data allow deeper analysis of the dust distribution. 



Cloud mass 

The cloud mass is derived from the extinction map using the following relation (Dickman 
1978): 

M=M)V^E^v'(^) (9) 

A.y 

I 

where a is the angular size of a pixel map, d the distance to the cloud, the mean molecular 
weight corrected for the helium abundance, and i represents a pixel map. With the dust- 
to-gas proposed by Savage & Mathis (1979), Nn/Ay = 1.87 x 10^^ cm^^ mag~^ {Nh = 
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Nhi + Nh^) and a distance of 580 pc, we obtain a mass of 4.5 x 10^ Mq. This value 
corresponds to about one fifth of the Orion cloud mass (Maddalena et al. 1986; Cambresy 
1999). Uncertainties come from the distance estimate for the cloud and the systematic error 
in the extinction. An uncertainty of a factor of two is believed to be a reasonable estimate. 
Figure 10 shows how the integrated mass distribution varies with the extinction within the 
cloud. Regions with visual extinction greater then 20 mag account for less than 3% of the 
total mass of the cloud. This suggests that the underestimation of the extinction in the 
cloud cores has only a limited consequence on the total mass estimate. 

Bally & Scovillc (1980) estimated the molecular cloud mass from CO observation to 
3-6 X lO'' M0 for a distance of 1 kpc. For 580 pc, it would correspond to 1-2 x 10^ M©. The 
discrepancy might come from the Nu/Ay value or from the conversion of CO density to H2 
density. 

3.3.2. Limitations of the method 

The cloud distance is the most critical parameter in the accurate application of these 
methods. The farther away a cloud, the lower is the spatial resolution and the larger the 
number of foreground stars. Here, we have mapped visual extinctions of 15 mag at a res- 
olution of 3' which corresponds to 0.5 pc {d — 580 pc). If the cloud were at 2 kpc, the 
3' resolution would correspond to 1.7 pc, equivalent to a study of nearby clouds such as 
Taurus, Chamaclcon or p Ophiuchus at a 40' angular resolution. This is the resohition of 
DIRBE/COBE and is probably close to the minimum acceptable to derive properties within 
a cloud. As the distance increases the number of foreground stars increases, as well, further 
degrading extinction maps. 

The 2MASS star density is about 4 times higher toward the galactic center than for the 
NAN direction. Clouds toward such very high density directions could therefore be mapped 
with a better resolution for a same distance. In contrast, the star density is about 3.5 
times lower toward the anti-center and dramatically decreases for higher galactic latitudes 
(density divided by 2 at 6 ~ 5°). Thus, the maximum distance for which molecular clouds 
can be mapped using 2MASS with this method cannot exceed 3 kpc toward very high density 
directions such as toward the galactic bulge. 



- 15 - 



3.3.3. Star Formation 

The technique described in section 2.3.2 allows one to identify star clusters by combining 
the color and the star count extinction maps. Since a cluster is an excess of stars compared to 
the local density, it should appear as a deficit of extinction in the star count map compared 
to the color map. To identify star clusters toward the NAN, we select areas for which 
Ay(color) — 74y(counts) > 6 ^ 1.7a and that contain two adjacent pixels with 74y(color) — 
Ay (counts) > 14 4.0(7. These numbers were chosen subjectively after examination of 
the full resolution 2MASS images. Nine cluster identifications are proposed in Table 1, and 
Figure 11 shows their position in the cloud. Figure 12 shows the 2MASS color images with 
the threshold contour overlaid. Among these nine clusters, only number 3b and 7 were 
already known (Herbig 1958; Wendker et al. 1983, respectively). 

Cluster 1 and 2. They are located on the edge of the central core of the dark cloud. 
The low density clustering of these red stars suggests either that the star formation occurs 
in the whole core (they are detected on the edge because the extinction is lower and the 
projection effect helps the detection at these locations) or that they are in the process of 
being dispersed. An upper limit for the star density excess can be obtained by assuming 
that the color map gives the true extinction and that the star counts underestimate the 
extinction because of the clustering. This yields only an upper limit because young stars 
are known to have an intrinsic infrared excess. For clusters 1 and 2 we find that the contour 
overlaid in Figure 12 corresponds to an excess of ~ 4 arcmin"^ and that the density reaches 
a maximum of ~ 20 arcmin^^ in their center. For comparison, clusters 5 and 6 have a similar 
value for the lowest contour but they both reach a maximum density 4 times larger of about 
80 arcmin"^. 

Clusters 3a and 3b. These two clusters are on the edge of the same globule. 3a contains 
faint red stars and 3b bright blue objects. Their relative position within the cloud may 
explain these colors (3a on the front side and 3b on the back side of the cloud). 3b was 
discovered by Herbig (1958). 

Cluster 4. It is located in a low extinction area and stars do not exhibit any infrared 
excess. It not clear whether this cluster is associated with the dark cloud; it might be a 
more distant open cluster. 
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Clusters 5 and 6. These are the two most obvious clusters from the star count and color 
comparison with Av'(color) — Ay(starcount) > 25 mag. They both lie on the edge of a highly 
extinguished globule. Cluster 6 is associated with an emission nebula that strongly supports 
the hypothesis of association with the cloud. 

Clusters 7 and 8. These are in the North America nebula (in the Hll region), in a low 
extinction area {Ay ~ 3 mag) and apparently not associated with any extinction core. The 
clustering is obvious for both clusters but there is no other proof of their association with the 
dark cloud at 580 pc. Cluster 7 has been observed by Wendker et al. (1983) who concluded 
it is located in the Perseus spiral arm, at a kinematic distance of 5-7 kpc, far behind the 
NAN. Cluster 8 is very compact compared to the others and that might indicate a larger 
distance. It contains 12 stars among which 3 are detected in JHKg, 3 in HKg and 6 in Kg 
only. 



4. Conclusion 

We have proposed a method to map extinction with 2MASS data using simultane- 
ously Kg star counts and H — Kg reddening. Finding a value of the extinction significantly 
larger from the count map than from the color map indicates significant foreground star 
contamination. Foreground stars can be removed individually in high extinction regions and 
statistically elsewhere. We stress the point that the foreground population dominates for 
high extinction and these regions would not be correctly analyzed if foreground stars are 
not carefully removed. On the other hand, a higher extinction value derived from colors 
than from counts indicates the presence of a star cluster. The adaptive method increases 
the angular resolution of the mapping around these objects. Color, especially H — Kg, is 
generally a better estimator of the extinction than star counts which have a higher statistical 
uncertainty. 

Towards the NAN, we estimate the foreground star density to 1900 ± 500 stars deg~^ 
and the visual extinction is mapped up to 35 mag. Angular resolution of the extinction map 
depends on the local stellar density and is 1.2', 3.7' and 7' for visual extinctions of about 0, 
20 and > 30 mag, respectively. Nine clusters are identified, only two of which were already 
known and another two of which correspond to a diffuse excess of red stars close the center 
of the cloud, suggesting an evolved young population or star forming activity in the whole 
core (~ 5 pc diameter). 

Near-infrared wavelengths are required to investigate such highly obscured fields. With 
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the help of the accurate 2MASS photometry, we have obtained results that supersede the 
previous work on the extinction in this area. The method developed here can be extended 
to other high extinction clouds with a significant fraction of foreground stars (i.e. located at 
more than 500 pc, in the galactic plane) that were not accessible before. The knowledge of 
the extinction toward the galactic plane is essential in order to constrain stellar population 
model parameters, such as the diffuse extinction and the giant star distribution, that are not 
well constrained at higher galactic latitudes (see appendix). The generalization of studies 
based on the comparison of the extinction (derived from star counts and color) with molecular 
lines and/or far-infrared and submillimeter emission is now possible with the 2MASS data 
releases. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

L. Cambresy acknowledges partial support from the Lavoisier grant of the French Min- 
istry of Foreign Affairs. 



A. Comparison with a galactic stellar distribution model 

The extinction calibration is based on the H — Kg color in the reference field. This 
field is assumed to be representative of the aggregate color over the field containing the 
NAN. Deviations from this color arc interpreted as extinction from the dust associated with 
the complex. To validate this procedure, we have compared the 2MASS star counts with 
the predicted counts based on galactic stellar distribution model, which predicts both the 
star counts and general interstellar medium extinction for the region of the Milky Way that 
includes the NAN. After incorporating the derived NAN extinction map with the model star 
counts, the extincted star counts should match the observed 2MASS star counts (§ A. 3). 

A.l. The Model 

The genesis of stellar distribution model comes from the Bahcall & Soneira (1980) optical 
star count model. Jarrett (1992) modified the Bahcall & Soneira model using the discrete 
formalism of Ehas (1978), Jones et al. (1981) and Garwood & Jones (1987), extending the 
model to the near- infrared (1-5 iim). The model includes the class III evolved giants, class 
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IV subdwarfs, class V main sequence, brown dwarfs (types L and T), and AGB populations. 
As with the Bahcall & Soneira model, the stars are distributed in two primary large-scale 
components: disk (exponential profile) and spheroid {R^^^ profile). Colors, luminosity, and 
number density per spectral type are based on the empirical data from Wielen (1974), Wielen 
et al. (1983), Koornneef (1983), Reid & Gilmore (1982), Hawkins & Bessell (1988), Bessell 
(1990), Bessell & Brett (1988). The interstellar extinction is applied as a smooth exponential 
function of galactic position, characterized by its own scale height and disk length. In total 
there are several adjustable parameters, including disk/spheroid scale lengths, scale heights, 
luminosity function and dispersion, and extinction/reddening laws. The model parameters 
were tuned using deep optical and infrared star counts (see below). 



A.2. Validation of the Model 

The optical portion of the model was tuned and validated using deep CCD observations, 
y < 24 mag, it! < 22 mag, of a diverse set of fields, encompassing both low and high stellar 
number density (Jarrett 1992; Jarrett et al. 1994). The near-infrared portion was tuned 
and validated using 2MASS star counts, for J < 15.5 and Kg < 14 mag, for fields ranging 
from 0° < |6| < 90° (Table 2). The typical field size was 25 deg^. The model performs well 
for |6| > 20°. An illustrative example of the performance for / = 90° b = 30° is given in 
Figure 13. At low galactic latitudes the model has small but significant differences compared 
with 2MASS, probably due to differential extinction from clouds along the line of sight (e.g., 
see \b\ < 5° in Table 2, corresponding to regions of the plane that are in close proximity to 
the NAN). A direct comparison with the NAN is given in § A. 3. 



A. 3. Comparing the model with the 2MASS NAN observations 

The galactic stellar distribution model adequately (~ 10-20%) predicts the number 
density and colors of stars located in the galactic plane, including areas near the NAN 
(Table 2). For the cloud itself, an additional source of extinction is required to match the 
model with 2MASS measurements. The histogram of the extinction we obtained using our 
extinction map allows us to redden statistically stars that are located more than 580 pc from 
the Sun (the nominal NAN cloud distance). We emphasize that this extinction component 
(attributed to the NAN cloud) is an addition to the nominal interstellar extinction component 
that is built into the model (§ A.l). 

Figure 14 compares the 2MASS data with the extincted-corrected model predictions. 
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The upper panels show the results for the cloud complex {I — 85°, b — —0.9°) and the 
lower panels show the results for the reference field (l — 86.2°, b = —2.9°). The lower left 
panel shows differential star counts for the reference region, where the solid and dashed 
lines correspond to the model and the diamonds denote the 2MASS measurements. Note 
that the evolved stars dominate the bright end {Ks < 11 mag) and disk dwarfs dominate 
the faint end. The H — Kg color histograms are shown in the lower right panel for the 
reference region. Both the star counts and colors suggest that the 2MASS measurements 
are well fitted by the model predictions (see also Table 2). The upper left panel shows 
the differential star counts for the region along the cloud hne-of-sight. The shape of the 
H — Kg color histograms are remarkably similar (right upper panel) , again suggesting that 
the model+cloud extinction adequately describes the stellar distribution toward the NAN. 
For comparison we show (upper right hand panel) the predicted stellar colors for stars located 
along the linc-of-sight that is directly opposite the NAN in the galactic plane (i.e., / = 85°, 
b = +0.9°; see Table 2). Here we have applied an ad hoc extinction. Ay = 1 mag at 4 kpc, to 
match the 2MASS star counts along this line of sight (Table 2, note (c)). This is consistent 
with the expected low extinction derived from CO emission along this line of sight with 
W{CO) = 11.08 K km s~^ (Dame et al. 2001) which corresponds to Ay = 1.1 mag. Star 
colors are relatively bluer and the total number of sources is nearly double to what is seen 
toward the NAN. 

The star count model suggests that the extinction in the reference region is probably 
quite small, consistent with HI observations and visual inspection of the POSS plates. More- 
over, the predicted number of foreground stars, 1900 ± 500 deg^, imphes a cloud distance, 
~ 485 pc, that is comparable (within the model uncertainties, ~ 100 pc) to the currently 
accepted distance to the NAN, ~ 580 pc. Given the overall high-quality consistency between 
the model predictions and the 2MASS measurements, we conclude that the reference field 
gives the correct zero point color in the whole NAN complex, and the derived extinction 
map is a fair representative of the dust column density attributed to the NAN. 
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Table 1. Coordinates Of The Star Clusters 



Id 


R.A. (J2000) 


Dec. (J2000) 


1 


20 56 02 


43 37 00 


2 


20 57 12 


43 48 30 


3a 


20 58 28 


43 56 24 


3b 


20 58 19 


43 53 37 


4 


20 56 13 


44 23 05 


5 


20 50 50 


44 25 00 


6 


20 53 42 


44 31 36 


7 


20 53 35 


44 47 20 


8 


20 54 14 


44 54 07 
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Table 2. 2MASS Fields Used To Constrain Model Parameters 





h(°) 




Performances 








star counts 


color 


90 


90 


Galactic North Pole 


nominal 


nominal 


90 


50 




nominal 


nominal 


90 


30 




nominal 


nominal 


90 


20 




nominal 


(a) 


270 


20 




nominal 


nominal 


90 


15 




nominal 


(a) 


90 


10 




(b) 


(a) 


180 


10 


Galactic Anticenter 


nominal 


nominal 


270 


10 




nominal 


nominal 


90 


5 




(b) 


(a) 


90 


2 




nominal 


(a) 


86 


-2.9 


reference field (Fig. 14, 
lower panels) 


nominal 


nominal 


85 


1 


opposite the NAN 


(c) 


(c) 



Note. — Star counts refer to counts at Kg — 14**^ mag; colors 
refer to the J — H, H — Kg and J — Kg mag differences. A 
nominal result means that the model star counts are within 
10% of the 2MASS star counts, (a): slight color excess in the 
2MASS J - Ks with respect to the model colors, (b): 2MASS 
i^s-band star counts in ~ 20% excess with respect to the model, 
(c): star counts and colors are nominal only if an additional (ad 
hoc) source of extinction is included in the model: Ay = 1.0 at 
4 kpc. 
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Ay (true) 

Fig. 1. — Effect of tlie foreground star contamination on the measured extinction versus tlie 
true extinction for different degree of contamination {X ~ jybackground^jy-totai^ different 
metliods {H — Kg color, Kg star counts and J star counts). 
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Fig. 2. — Histogram distribution of the 2MASS magnitudes in the NAN (J: dashed line, H: 
solid hne and Kg. dashed-dotted hne). Slopes of the luminosity functions are 0.36, 0.34 and 
0.34 for J, if and Kg., respectively. 




Fig. 3. — Digitized Sky Survey optical image (left) and 2MASS J — Kg map (right) of the 
NAN. The square shows the area studied in the paper. The reference field is used to set the 
zero point of extinction. Coordinates are expressed in the equatorial J2000 system. 
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Fig. 4. — Color color diagram for the reference field (top) and the NAN (bottom) directions. 
RL85 and CCM89 are the reddening vectors for Rieke & Lebofsky (1985) and Cardelli et al. 
(1989) extinction laws, respectively. Dwarf and giant sequences are overlaid (Bessell & Brett 
1988). 
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Fig. 5. — Velocity components toward the NAN (solid line) and the reference field (dashed 
line) from the Leiden/Dwingeloo survey. The vertical dotted lines show the velocity range 
of the NAN (in the Local Arm). 




Fig. 6. — Extinction maps from Kg star counts (top) and H — Kg color (bottom) for a 
64' X 24' field centered on the most obscured region of the dark cloud (q;j2ooo = 20^56™51*^, 
(5j2ooo = 43°42™33^). Black is for high extinction and white for low extinction. Extinction 
from star counts reaches 30 mag whereas the extinction from reddening at the same position 
is mag due to the foreground star contamination (see § 2.3.2). 
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Fig. 7. — Distribution of foreground stars with smoothed extinction contours overlaid (con- 
tour values are 5, 10 and 15 Ay). Extinction is derived from H — Kg color. 
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Fig. 8. — Extinction maps derived from i^., star counts (left) and from H — Ks color (right) 
corrected for foreground stars. Coordinates are expressed in the equatorial J2000 system. 
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Fig. 9. — Extinction map of the NAN obtained by combination of the H — Kg color map 
with the Kg star count map. 
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Fig. 10. — Cumulative enclosed mass in iso-extinction contours versus visual extinction. 
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Fig. 11. — Position of the nine clusters (contours correspond to Ay (color) — Ay (counts) > 
6 mag) represented on the combined extinction map of the cloud. 




Fig. 12. — JHKg color images of the clusters selected by comparison of the local stellar 
density with the median reddening. The width of each field is specified in the lower left 
corner. Contour for Ay (color) — Ay (counts) > 6 mag is overlaid when the clustering is not 
obvious. 
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Fig. 13.— Star counts and colors for / = 90°, b = 30°. A 5° x 5° 2MASS field, centered on 
galactic position / = 90° b = +30°, is depicted with A'^-band star counts and J — H, H — Ks 
and J — Kg colors. The color histograms are restricted to 6 < < 14 mag. Overlaid are 
the model predictions. 
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Fig. 14.— X^-band star counts and H - Ks colors for the NAN (/ = 85°, b = -0.9°), 
upper panels, and the reference region (/ = 86.2°, b = —2.9°), lower panels. The 2MASS 
measurements are depicted with diamond symbols. Overlaid (connected lines) are the model 
predictions. For comparison, the upper right panel includes the predicted H — Kg colors for 
stars located along the line-of-sight that is directly opposite the NAN in the galactic plane 
{I — 85°, h — +0.9°; see also Table 2). NAN model includes extinction map. 



